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We describe a method to probe the spectral fluctuations of a transition over broad ranges of 
frequencies and timescales with the high spectral resolution of Fourier spectroscopy, and a temporal 
resolution as high as the excited state lifetime, even in the limit of very low photocounting rates. 
The method derives from a simple relation between the fluorescence spectral dynamics of a single 
radiating dipole and its fluorescence intensity correlations at the outputs of a continuously scanning 
Michelson interferometer. These findings define an approach to investigate the fast fiuorescence 
spectral dynamics of single molecules and other faint light sources beyond the time-resolution 
capabilities of standard spectroscopy experiments. 



INTRODUCTION 

Chromophores embedded in a condensed medium inevitably exhibit time-dependent, fluctuating optical properties 
reporting on the dynamics of their nanoscale environment. Studied on ensembles of molecules [l|, these fluctuations 
illuminate the complex dynamics of a broad range of disordered host systems such as low-temperature glasses, proteins 
and liquids, over timescales extending from femtoseconds to hours (2,y|- Observed at the single emitter level, chro- 
mophore fluctuations reveal surprisingly varied and complex dynamical phenomena kept hidden in ensemble-averaged 
experiments jj|, whose understanding is crucial for developing applications such as single molecule probes in biophysics 
or single photon sources in quantum information processing pj. 

Single molecule spectroscopy reaches its fullest potential when combined with high time resolution, so as to resolve 
the fast processes and temporal heterogeneities of any given isolated emitter. Obtained down to timescales shorter than 
the excited state lifetime for single emitter intensity fluctuations ja|, high time resolution proves much more difficult 
to achieve when probing spectral fluctuations, due to the flnite integration time (typically larger than milliseconds) 
necessary to collect enough photons to measure the spectrum of the emitter |2| . 

The first attempt to achieve fast single molecule spectroscopy was undertaken by Plakhotnik, who demonstrated 
that time resolution could be improved by averaging the autocorrelation functions of many individual fast scan spectra 
provided by single molecule laser spectroscopy (a,|2l|- The time resolution of this intensity-time-frequency correlation 
(ITFC) method is however technically limited to the order of milliseconds by the finite scan rate of the laser, and 
no application followed this pionneering work. As a result, our current understanding of fast spectral fiuctuations in 
single emitters is so far mostly inferred from inhomogenously broadened ensemble experiments. 



/^ 



single emitter 
fluorescence 



" u 



la(t) 



vt 



>'~ 



il!) start 4 
Pw^\Jl stop o- 



V 



lb(t) 



beam splitter 



Fig. 1 - Single molecule photon-correlation Fourier spectroscopy setup. Starting from an initial optical path 

difference Si, the intensity correlation function g{T) of the output intensities Ia{t) and Ih{t) is measured during a 

continuous scan of the interferometer at a velocity V. Repeating this procedure for various values of Si provides the 

time-resolved frequency fiuctuation spectrum PriC) of the emitter. 



In this article, we describe an approach overcoming these limitations to explore the spectral dynamics of a single 
transition at both high temporal and high spectral resolutions with standard photocounting equipment, and without 
requiring mode-locked laser sources for excitation as in ITFC. The method is based on the observation that spectral 
fluctuations under continuous excitation can be directly encoded in the intensity correlation functions measured at the 
output of a continuously scanning Michelson interferometer by a pair of photodiodes followed by a photon correlation 
counting board [Fig. ^. The setup is investigated theoretically, simulated numerically, and compared to standard 
spectroscopy experiments. 

THEORETICAL BACKGROUND 

Intensity correlation functions at the output of a scanning interferometer 

The basis of our approach is the existence of an intimate relation between the spectral fluctuations of a transition and 
the intensity correlation function of its radiated field after transmission through a scanning Michelson interferometer. 
We first derive this relation in the case of a monochromatic transition undergoing random stationnary fiuctuations 
Suj{t) around a frequency uiq — {uj{t)) over a mean-squared range a^ — {SLu{t — 0)^), and later generalize to fluctuating 
transitions with any lineshape. 

Starting from an initial optical path difference Si at time t = 0, the mirror travelling at a constant velocity V creates 
a time dependent optical path difference 6{t) = Si + 2Vt between the two arms of the interferometer. After entering the 
interferometer, frequency fluctuations propagate in the arms and recombine on the beamsplitter. This recombination 
can be described as instantaneous, as the retardation effects (of order Si/c) induced by optical path differences Si 
shorter than the maximum possible coherence length cTi of the emitter remain neglible over the timescales r > Ti 
considered throughout this paper (i.e. 5i/c <^ t. See also (ill). Hence, we have : 

hit) ex l + cos[{2Vt + S,)uj{t)/c] 

h{t) ex l~cos[{2Vt + S,)uj{t)/c]. ^ ' 

These oscillating intensities feed a photon counting board integrating photons between i = up to a time t = T for 
computing the time-averaged intensity correlation function between the two outputs of the interferometer 



9{r) . ^(!M±lL, (2) 

h{t) h{t + T) 

where ... denotes time-averaging from t = up to i = T. The integration time T is chosen so that the corresponding 
change in optical path difference A = 5{T) — (5(0) = 2VT spans a large number of fringes, reducing the time-averages 
in Equation 12 to : 

5(t) ^1-^1 cob{2looVt/c + a{t)5^/c)dt (3) 

with a{t) — 26Lo{t + T)VT/5i + (1 + 2Vt/6i)CT{t), where Cr(0 — ^^{^ + ■^) ~ 5oj{t) denotes the random frequency 
fluctuation observed between times t and t + r. 

This general expression of g{T) simplifles with an appropriate choice of the scanning parameters V and 5i, since the 
flrst term oi a{t) becomes negligible in g{T) at slow scanning velocities, when 25ijj{t-\-T)VT jSi ~ 2(jVT/5i <C c/i5i, i.e. 
as the change in the optical path 2Vt occuring over the timescale r under investigation remains small compared to 
the coherence length A = c/a of the emitter. The quantity a{t) then reduces to a{t) = C,r{t) in Eq. Elwhen 5i :s> A, 
and the intensity correlation function becomes : 

5(t) = 1 - ^ / cos(2wo't^T/c + Cr{t)5,/c)dt. 

For integration times T much larger than the typical timescale over which spectral fluctuations occur, time averages 
can be replaced by the ensemble-average over the distribution Pt(C) of all the possible realizations of the random 
variable C = Ct(^ = 0), so g{T) now becomes : 

9{r) = ^-\j cos{2ujoVt/c + Q5^/c)pr{C,)dC, 



which can be recast in a more compact form as : 

5(t) = l-\ cos(2u;oVT/c)¥T[pr{0]s^/c (4) 

when the fluctuation process is time-reversal invariant, which imposes Pr(C) =Pr(~C)- 

Hence, provided the two conditions i) 2Vt ^ c/a and ii) 5i^ IS. are fulfllled, the time-averaged intensity correlation 
function g{T) measured at the output of a scanning Michelson interferometer oscillates with the frequency 2loqV/c 
at which fringes oscillate on the photodetectors, and with an amplitude given by the value of the Fourier transform 
¥T[pr{C,)]{0) at 9 — 5i/c, where Pr(C) is the distribution of frequency shifts. 

General case : narrow transition with arbitrary lineshape 

The result obtained above holds for fluctuating transition of any (narrow) lineshape, integration over the flnite 
linewidth simply changing the distribution of frequency shifts Pt{C) into the more general expression : 

^+00 

PriO = { Stiuj)st+ri^ + Odu;}, (5) 



where (...) denotes ensemble averaging over all possible realizations of spectral fluctuations, and St(w) is the time- 
resolved emission spectrum of the transition |:2| : 



1 /"-t-oo , 

where Ti is the excited state lifetime of the transition. Remarkably, Eq.^identifles with the time-frequency spectrum 
provided by past ITFC experiments, and is the Fourier transform (in () of the echo obtained in three-pulse photon 
echo spectroscopy |i2, I^J . 

Photon-correlation Fourier spectroscopy 

The key point is that at timescales t shorter than the periodicity ttc/ujoV at which the fringes oscillate on the 
photodetectors, condition i) (cf. 2.1) is automatically fulfllled for any narrow spectral fluctuations (i.e. when a ^ ujq) 
and Eq. 0] can be reverted as : 

p,(C) = 2FT-i[l-g(r)]^=2W5.- (7) 

The distribution of the spectral fluctuations PriC) oi the emitter can therefore be determined directly from the 
measurement of its fluorescence intensity correlation function g(r) at various optical path differences 6i of the scanning 
interferometer. 

Determining the distribution Pr(C) of an emitter rather than its fluorescence spectrum Sf(w) suggests an original 
approach to investigate its fast spectral dynamics, very much in the same way that measuring an intensity autocor- 
relation function instead of an intensity time trace transformed the study of photon statistics on short timescales in 
faint light sources such as distant stars or nanoscale emitters [ll|. This can be understood by taking a closer look at 
some of the general properties of the distribution PriO- When r — > oo, the spectra St(w) and St+ri^) are statistically 
independent, and so Pt(C) reduces to the autocorrelation of the time-averaged (inhomogenous) linewidth s{uj — ujq) of 
the transition accessed by any standard fluorescence spectroscopy experiment. On the contrary, in the limit r ^ 0, 
PriC) is given by the autocorrelation of the time-resolved (homogenous) Lorentzian lineshape of purely radiative width 
Tj~^, as no fluctuation has time to occur within delays r = 0. The distribution Pt(C) hence naturally bridges the gap 
between single molecule fluorescence trajectory analysis and ultra-fast spectroscopy ensemble experiments. 

The time-frequency distribution Pt{C) also has the property of being time-independent for any stationnary dis- 
tribution of fluctuations. The measurement oi Pr{C) can therefore be made by integrating photon coincidences over 
long durations to improve signal-to-noise ratio without degrading temporal and spectral resolutions. Derived within 
the frame of classical electrodynamics, our analysis extends down to timescales Tyain as short as the excited state 
lifetime of the emitter Ti , readily accessible by the modifled Hanbury-Brown Twiss photon correlation detection setup 
described Fig. 1 [li|. This approach - which we call photon correlation Fourier spectroscopy (PCFS) - therefore com- 
bines the high spectral resolution Cmin ~ c/5i of Fourier spectroscopy with a high temporal resolution down to the 
Fourier-transform limit TniinCmin '^ 1, beyond the current capabilities of standard single molecule experiments. 
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Fig. 2 - Photon correlation spectroscopy of a single static (left) or switching doublet (right), (a) Intensity 

correlation function at various delays 6i . The scatter plots are numerical simulations for an emitter detected with an 

intensity / = 50 kHz. (b) Evolution oi gij) with 6i for r=4 ns (o), 2.5 /xs (D), 10 /is (o), 160 /j.s (•), depending of the 

optical delay bi where the measurement was performed, (c) Corresponding fluctuation distribution Pt(C) (°)- (d,e,f) 

Same as in (a,b,c) for the switching doublet. Solid lines are the theoretical expectations corresponding to the 

simulation parameters (see Table 1). 



NUMERICAL SIMULATIONS AND DISCUSSION 

Numerical simulations were performed to investigate the validity of PCFS for exploring fast spectral fluctuations 
and highlight some of its instrumental properties. All simulations were made assuming an excited state lifetime 
Ti = 1 ns, and fast random spectral fluctuations over a total range (SA — 2AoO'/cJo — 1 nm around an average 
wavelength Ao = 600 nm, with an exponential frequency correlation function Cij') = {6uj{t)5ijj{t + r)) of correlation 
time Tc = /q C{t)dt/a'^ = 5 fis. The choice of these numerical parameters corresponds to typical values as encountered 
in standard single molecule spectroscopy experiments for emitters such as molecules or semiconductor quantum dots. 

PCFS measurements were numerically simulated with a speed of the translation stage set to y = 30 /im/s to 
fulfill condition i) {V <C 1.8 mm/s) over timescales r < 100 ms. The photon arrival times at the entrance of the 
interferometer were drawn according to the Poissonian statistics of a light beam of intensity / = 5 x 10'* photons/s. 
This corresponds to a situation where spectral fiuctuations occur at timescales Tc — 5 fj,s much faster than the average 
delay I~^ =20 /is between two successive photodetection events, and so are completely averaged out in standard single 
molecule spectroscopy experiments. The photons were statistically directed towards either photodiode depending on 
their wavelength and arrival time in the interferometer. Intensity correlation functions were calculated for various Si 
by integrating photons over 30 fringes over 500 scans, corresponding to a total acquisition time of 5 min per intensity 
correlation function. 



PCFS and spectral fluctuation dynamics analysis 



We first investigated the outcome of a PCFS experiment for a transition undergoing discrete spectral fiuctuations, 
with a frequency uj{t) switching between two values ui and 102 = uji + fl (here corresponding to wavelength jumps of 
±1 nm) as a random telegraph signal S{t) [Fig. 2, Tab. 1] as encountered in the study of chromophores interacting 
with the two- level systems in glasses at low temperature [3|- The transition frequency fiuctuations were generated 
with exponentially distributed waiting times of rate ki and fe with fci = /c2 (and Tc — [fci + fc2]~*). Observed in 
standard spectroscopy, the transition spectrum would appear as predicted by the Anderson-Kubo lineshape theory 
[131, i.e. a doublet of separation fl = 2a centered in [cui + W2]/2, indiscernible from the spectrum of a static doublet 
of transitions at frequencies wi and a'2 with similar intensities. 

Figure 2 shows the result of PCFS experiments on a static doublet (left), and on the switching transition (right). 



5co(t) 
Sudden jumps Pi^ 

-P2Q 



2 



Q 



Uj{t) = CJO + SlS(t) 

C{T) = Q.'^{S{t)S{t + T)) 



g{T) = 1-5 cos(2LJol^r/c)e-'''/"^i \Ao + An cos(™,/c)] 
' Ao = (piP2)'b? +pi + 2C(r)/n2] 



where , 

An = 2(pip2)'[C(0)-C(r)]/n2 



5(0(t)i 
Gaussian fluctuations 2ot ^"V^yArN ,VV» t 



Prob(<5aj) oc e-*" ^^.r 
C(r) = cr'e-"/"^ 



g{r) = 1 - i cos(2a;ol^r/c)e-''/'=^iF,((5./c) 
where Fr(t) = e-io(«-«')l2C{t')-C(t'+r)-C(t'-T)]dt' 



Tab. I - Theoretical expression of the intensity correlation function ^(t) measured in PCFS for discrete and 

continuous spectral fluctuations. C{t) = {Suj{t)SLu{t + r)) is the frequency correlation function of the fluctuations. 

Pi=i,2 denote the fraction of time spent by the transition in states 1 and 2 respectively. 

As seen in Fig. 2(a)(d), in both cases, the shape of the intensity correlation function strongly depends of the delay Si 
where it was measured. Repeating the measurement of g(r) over different optical path differences Si, we determined 
the evolution of ^(r) with Si [Fig. 2(b)(e)], from which the distribution Pt{C) was extracted with Eq.[3[Fig. 2(c)(f)]. 
For a static spectrum, photons do not exhibit spectral correlation, and so PriC) reduces to the autocorrelation of 
the time-averaged spectrum measured in standard spectroscopy (here a doublet), i.e. a triplet of lines of intensities 
{1/4,1/2,1/4} at frequencies {— fi, 0,+r2} independently of the timescale r, as observed in the experimental results 
shown Fig. 2(c). 

For the switching transition, this pattern is only preserved over timescales where fluctuations are uncorrelated, i.e. 
when r > Tc = 5 /is [Fig. 2(f)], and breaks down as soon as r '^ Tc, when the sidebands oi Pr{C) in C = ±fl decay 
progressively as r decreases, asymptotically leaving us with the autocorrelation of the time resolved spectrum of the 
transition as r — > - here a Lorentzian of width Tj~^. Interestingly, the calculation of g(r) from Eq. 14161 (cf. Tab. 1) 
indicates that the correlation function of the random telegraph signal governing the frequency switching process is 
directly encoded in the decay of the sidebands amplitude Aq with r. 

The insets in Fig. 2(c)(f) - showing the correlation functions C{t) extracted from the decay of the measured sideband 
amplitudes An - conflrm this prediction. Indeed, they indicate respectively that the static and dynamic doublet exhibit 
null and exponentially decaying correlation function, as expected. Illustrated here on binary spectral jumps, the ability 
of PCFS to investigate spectral fluctuations dynamics down to the excited state lifetime of the emitter can also be 
exploited on a transition coupled to a collection of flipping two- level systems, for example to determine the detailed 
physical properties of the latter (energy splittings Jl, switching timescales Tc, etc.) [2|. 

High-resolution spectroscopy beyond temporal inhomogenous broadening 

PCFS also opens a perspective for performing high resolution spectroscopy despite the presence of broad, fast, 
continuous spectral fluctuations, as reported on most emitters - often down to low temperatures (4.2 K) |lj]. This 
point is illustrated Fig. 3, showing the result of a PCFS experiment simulated for a doublet undergoing fast stationnary 
Gaussian fluctuations over a range a much broader than the doublet separation fl (here a — 5il) , with an exponential 
correlation function C{t) (Ornstein-Uhlenbeck fluctuation process) with Tc = 5 /xs. 

Observed in conventional spectroscopy, the transition would appear as a Gaussian lineshape of width a, and the 
existence of the underlying doublet would remain unnoticed. PCFS, in comparison, reveals a completely different 
pattern. Here again, the intensity correlation functions were found strongly dependent oi Si [Fig. 3(a)], and, measured 
for different values of Si, provided the distribution Pt{C) [Fig- 3(b)(c)]. As seen in Fig. 3(c), the distribution PriC) - 
broad (here a Gaussian of FWHM=2cr) over durations t > Tc - progressively narrows when t < Tc. This is consistent 
with the fact that for photons separated by delays t shorter than the fluctuation correlation time Tc, fluctuations are 
increasingly seen as "frozen" as t ^ 0. 

Due to this fluorescence line-narrowing effect, the doublet can be resolved over short timescales, as seen from the 
quasi-periodic oscillations (of periodicity 2Trc/fl in Si) appearing in the intensity correlation function gir) when r ~ 5 
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Fig. 3 - Photon correlation spectroscopy of a doublet of separation fi undergoing Gaussian stationnary fluctuations 

of correlation time Tc = 5^j,s, over a spectral range a = 5fl (corresponding to SX = Inm). (a) Intensity correlation 

function at various delays Si obtained from numerical simulations when the emitter is detected with an intensity 

/ = 50 kHz. (b) Evolution of g(r) with Si for r=2 ns (D), 40 ns (x), 640 ns (o) as observed from the measurement of 

gir), depending of the optical delay Si where the measurement was made. At short timescales (r < 10 ns), 

oscillations of periodicity 2nc/Q are observed, as the doublet becomes resolved, (c) Corresponding fluctuation 

distribution PriO- ^t short timescales, a triplet appear, i.e. the doublet is resolved. Solid lines are the theoretical 

expectations corresponding to the simulation parameters (see Table 1). 



ns [Fig. 3(b)], which translate into the lineshape of pr<5ns(C) expected for a static doublet, i.e. a triplet of intensities 
(1/4,1/2,1/4) at frequencies (— ri, 0, fl) [Fig. 3(c)]. A corollary of the fluorescence Hne-narrowing effect is an increased 
coherence length of the emitter, coherence in the transition emission (i.e. ^(t) < 1 in [Fig. 3(b)]) remaining visible 
over short timescales t ^Tc although we have 5i> b mm ^ A, i.e. even if the optical path difference is much larger 
than the bare coherence length A — 120 /im of the emitter, when no fringe would be observed in standard Fourier 
spectroscopy. 



EXPERIMENTAL PROPERTIES 



From a technical standpoint, the intensity correlation measurements at the basis of PCFS can be considered as 
an intensity homodyne detection, where the oscillating intensity detected by a photodiode is demodulated by the 
oscillating intensity detected on the other photodiode. A first consequence of this observation is that PCFS directly 
provides the envelope of the Fourier transform interferogram of the radiation field, without the complex demodulation 
schemes usually involved in scanning Fourier spectroscopy. 

Secondly, the self-demodulation process implies that the shape of g(r) is robust against fluctuations in the scanning 
velocity V (as caused by stick-slip and vibrations in the translation stage) , and is independent of the exact average 
frequency ojq of the transition, making PCFS intrinsically insensitive to rare, large spectral jumps which often limit the 
measurement time of laser spectroscopy experiments (e.g. the PCFS inversion formula (Eq. 8) is indeed independent 
of V and ojq). Numerical simulations (not shown) confirmed this analysis. For example, random velocity fiuctuations 
of 30% were found to have no significant impact on any of the simulated results presented in Fig. 2 and Fig. 3. 

Finally, we note that PCFS also offers high time and high spectral resolution over a broad range of frequencies and 
timescales, contrasting with laser spectroscopy, which - because of its very scanning nature - only provides high time 



and spectral resolutions simultaneously in the limit of vanishing spectral ranges. 

CONCLUSION 

Replacing the beamsplitter of a Hanbury-Brown Twiss detection system by a scanning Michelson interferometer 
allows the measurement of spectral fluctuations of a transition at high spectral resolution, down to timescales as short 
as the transition excited state lifetime, which opens unexplored possibilities for studying the fast fluorescence spectral 
dynamics of single dipoles in a range of contexts. 

In solid state physics, this photon-correlation spectroscopy method suggests an approach to investigate the relation 
between lineshape broadening and decoherence in molecules and semiconductor quantum dots, as well as the dynamical 
interactions of these systems with optical fields and their nanoscale environment. Implemented in a fiuorescence 
correlation spectroscopy (FCS) experiment, PCFS might also provide some insight into the spectral dynamics of 
nanoscale emitters under the influence of chemical reactions, conformational changes or intermolecular interactions in 
liquid environments. 
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